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ABSTRACT 

The low-luminosity black hole Sgr A* provides a testbed for models of Radiatively Inefficient Accretion 
Flows (RIAFs). Recent sub-millimeter linear polarization measurements of Sgr A* have provided evidence 
that the electrons in the accretion flow are relativistic over a large range of radii. Here, we show that these 
high temperatures result in elliptical plasma normal modes. Thus, polarized millimeter and sub-millimeter 
radiation emitted within RIAFs will undergo generalized Faraday rotation, a cyclic conversion between linear 
and circular polarization. This effect will not depolarize the radiation even if the rotation measure is extremely 
high. Rather, the beam will take on the linear and circular polarization properties of the plasma normal modes. 
As a result, polarization measurements of Sgr A* in this frequency regime will constrain the temperature, 
density and magnetic profiles of RIAF models. 

Subject headings: accretion, accretion disks — black hole physics — Galaxy: center — radiation mechanisms: 
nonthermal — radiation mechanisms: thermal 
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1. INTRODUCTION 

Radiatively Inefficient Accretion Flows (RIAFs) are models 
of black hole accretion at rates much smaller than the Edding- 
ton rate. In this regime, numerous physical effects reduce the 
ability of the infalling gas to radiate, thereby robbing accre- 
tion of its capacity to convert gravitational energy into light. 
For example, in an early RIAF model known as an Advection 
Dominated Accretion Flow (ADAF) the gas density is low 
enough that the ele ctrons and ions thermally decouple (e.g., 
IShapiro et al.ll 19761) . so that the accretion energy is carried by 
the ions into the black hole with little escaping as radiation 
dIchimarull977HRees et aUflOSllNaravan & Yil lT995). Later 
work found that these low density flows resulted in gas that 
is unstable to large radial circulations (Ouataert & Gruzinoy 
UpOOa; Igumenshchev et al. 2000) and only mar ginally bound 
to the black hole tBlandford & Begelmanlll999l) . Thus, given 
a tiny perturbation, only a fraction of the gas that is fed to 
the outer regions of the accretion flow will actually success- 
fully fall into the black hole, further reducing the radiative 
efficiency. 

The benchmark against which all RIAF models are tested 
is the black hole at the center of the Galaxy, Sgr A*, which 
radiates at 3 x 10^^ of its Eddington luminosity, and has a 
well determined spectral energy distribution from radio wave- 
lengths to X-ray energi e s in both quiescent and flaring states 
(iMelia & Falckd IIOOU iBaganoff et aP llOOl iGenzel et alJ 
l2003h . Seve ral different RIAF models fit the spectrum of 
Sgr A* (e.g.,'Nai-avan et al." 1995': 'Ouat aert & Naravanl l 19991: 
lOzel e t al. 2000; Yuan et al. 2003), resulting in different solu- 
tions to the desired physical parameters such as the run of tem- 
perature, magnetic field stre ngth, and den s ity wi th radius. The 
constraints tightened when lAitken et al.l ( |2000|) reported the 
tentative detection of linear polarization in the sub-millimeter 
spectrum of Sgr A*. Agol (2000) and Ouataert & Gruzinov 
( I20b0b) then showed that the higher density ADAF solu- 
tions would likely depolarize any emission in this wavelength 
range, and speculated that only the hotter, more tenuous mod- 
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els are viable candidates. Now, with several successful ob- 
servations, linear polarization measurements are being used 
to constrain the a ccretion rate a nd density profile of the flow 
onto Sgr A* (e.g.. lMarrone et a l. 2006, 2007). 

This paper discusses how the relativistic plasma in RIAF 
solutions will effect the polarization properties of any propa- 
gating radiation. In particular, we show that observations of 
both linear and circular polarization at multiple frequencies 
can provide important constraints on many of the key predic- 
tions of RIAFs. 

2. ROTATION AND CONVERSION OF POLARIZED LIGHT 

Radiation propagating through a plasma pierced by a uni- 
form external magnetic field will travel via two wave modes 
that have orthogonal polarizations and different phase veloc- 
ities. The polarization properties of the radiation are thus al- 
tered as it passes through the medium. The best known exam- 
ple of this effect is Faraday rotation, where a linearly polar- 
ized wave has its plane of polarization rotated during passage 
through a cold (i.e., 9^1, where 9 = kTe/mc^ and Tg and 
m are the electron temperature and mass) plasma with circu- 
larly polarized normal modes. In general, however, both the 
plasma modes and radiation have elliptical polarizations giv- 
ing rise to generalized Faraday rotation. 

An effective way to visualize the result of this general 
case is to construct a three-dimensional space using the 
Stokes parameters as the Cartesian coordinates {Q, U, V)/I. 
The unit sphere in this s p ace is the Poincare spher e (e g. , 
iKennett & Melrosd 119981: iRuszkowski & BegelmanI l20()2h 
and is shown in Figure [T] The surface of this sphere de- 
notes all possible polarizations of the radiation and plasma 
wave modes. The poles and the equator of the sphere indi- 
cate either pure circular or linear polarization, respectively. 
The two orthogonal normal modes of a plasma are found 
as a line passing through the center of the sphere and one 
mode is drawn as the vector Cl in Figure [T] A polarized light 
ray N passing through this medium will rotate around the 
wave mode at a constant angle (shown in Fig. [T] as the circle 
on the sphere's surface) changing its polarization properties. 
The amount of this generalized Faraday rotation depends on 
the relative orientation between the polarization vector of the 
propagating radiation and of the normal modes of the inter- 
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Fig. 1 . — The Poincare sphere is the unit sphere in a space defined by the 
three orthogonal Stokes parameters normalized by the total intensity /. Any 
polarized radiation can be described by a point on this sphere with the poles 
and equator denoting pure circular or linear polarization, respectively. The 
vector Q shows one of the normal modes of a magnetized plasma. Polarized 
light, denoted by JV, will rotate around f2 at a constant angle as it passes 
through the plasma. Classic Faraday rotation con'esponds to the case when 
the modes are coincident with the poles. 

vening medium. The precession of the polarization vector can 
be analyzed by decomposing it into one component, iVn, that 
is parallel to £7 and one component, iVx that is perpendicular 
to it. The parallel component is unaffected by the propaga- 
tion through the medium, whereas the perpendicular compo- 
nent undergoes rapid rotation on the plane perpendicular to 
51, which results in a cyclic conversion between linear and 
circular polarization (Pacholczvk 1973). This mechanism for 
generating circular polarization has long been considered as a 
possible origin f or the observed circular polarizat ion from ra- 
dio quasars (e.g.. lRuszkowski & Begelmanll2002l) . It also has 
been applied to Sgr A* i n the context of a jet emission model 
dBeckert & Falckdl2002h . 

In a cold plasma the wave modes are coincident with the 
poles, so polarized radiation will rotate around a line of con- 
stant latitude on the sphere resulting in the classical Faraday 
rotation. If such a cold plasma produces a large number of 
rotations around the normal modes (i.e., it has a large rota- 
tion measure), then the radiation may be linearly depolarized 
as the linear polarized vectors are averaged out over different 
lines of sight. Interestingly, this does not occur for general- 
ized Faraday rotation. Even if the perpendicular component of 
the polarization vector suffers very large rotations around 51, 
the emerging radiation will still have a net elliptical polariza- 
tion described by the component TVn of the incident radiation 
parallel to the normal mode. Therefore, there will always re- 
main some linearly polarized component to the radiation and 
it is not generally true that a large rotation measure will lead 
to a linearly depolarized signal. 

To begin to identify the polarization signatures of a RIAF 
we first must calculate the ellipticity of the normal modes in 
the accretion flow. The axial ratio of the elliptical wave modes 
in a plasma can be written as (e.g., Melrose 1997); 

(a" _ ^ [(all _ ^22)2 ^ 4^12^211 1/2 
T± = - ^ „. 1, (1) 

where e*-' = + {^-Kcj uP'^a^^ is the plasma dielectric ten- 



Lo = 2ttv, where v is the frequency of the propagating radia- 
tion. The general expression for e'^ for a thermal magnetized 
plasma of any temperature was derived by Trubnikov ( 1958), 
and can be expressed in an analytic form in either the very 
low temperature {9 <C 1) or very high temperature {9 ^ 1) 
limit. If T = or cx), then the wave modes are linear and 
lie along the equator of the Poincare sphere, and significant 
conversion between linear and circular polarization will oc- 
cur. If T = ±1 then the modes are circular and only standard 
Faraday rotation is possible. 

Examination of equation [T] shows that the key parameter 
in determining the axial ratio of the linear modes is (aH — 
q;^^)/q!^^. For a plasma that has both a cold component with 
number density ricoid and a relativist ic {9 ^ 1) component 
with density n^^i then dMelrosdl 19971) 
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where is the angle between the magnetic field and the di- 
rection of propagation, and wb = eB/mc is the cyclotron 
frequency of an electron in a magnetic field of strength B. 
These two expressions are valid only for lo ':s> lo-q- For a 
non-relativistic plasma, the wave modes are independent of 
temperature, so the 9 in equations |2] and [3] refers to the tem- 
perature of the relativistic plasma. 

Consider now a thermal plasma with a temperature kTe 
that can range from non-relativistic to ultra-relativistic, such 
as within an RIAF (e.g.. Yuan et al. 2003). We can use the 
above equations to determine how the wave modes of this 
plasma will vary with temperature and thus its effects on po- 
larized radiation. Without loss of generality, the distribution 
of energies in this plasma can be written using the relativistic 
Maxwellian: 



(kT^ / •mc'^)K2{mc^ / kTe) ' 



(4) 



where 7 is the Lorentz factor, j3 is the electron velocity in units 
of c, and K2 is the modified Bessel function. For any value 
of the electron temperature, a critical Lorentz factor 7crit can 
be defined, above which a fraction of this plasma will always 
be considered ultra-relativistic. Integrating equation |4] above 
7crit allows us to split the thermal plasma into a relativistic 
component with density rirci and a cold component with den- 
sity TT-coid = 1 — 'Jrci- Re-defining 9 = (E) /mc^, where (E) 
is the mean energy of the plasma above 7ci it, the three quanti- 
ties (a^^ — a^^), a^^, and T can all be calculated as functions 
of the electron temperature. 

Figure|2]shows the result of such a calculation and plots the 
axial ratio of the normal modes T as a function of plasma tem- 
perature for an observed frequency of 100 GHz and a mag- 
netic field of 100 G. A value of 7ciit = 10 was used for this 
and all subsequent calculations. While approximations have 
been made to obtain this result, we see the expected behavior 
that the modes are purely circular (T = —1) for tempera- 
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It is in this regime that ADAFs reside (e.g.. 



lOzel et al]l2000l) and so they will only cause classical Fara- 
day rotation. The solid line in Figure |2] shows the value of T 
assuming a magnetic field that is tilted at <j) — 45° to the di- 
rection of propagation, while the dashed line shows the results 
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Fig. 2. — Axial ratio of the polarization ellipse (T) defined by the normal 
modes of a magnetized thermal plasma with temperature fcTe . A ratio of — 1 
indicates modes which are purely circular (i.e., they lie along the poles of 
the Poincare sphere), while a value of denotes modes which are entirely 
linear (i.e., they lie along the equator of the Poincare sphere). The plot shows 
the value of T assuming an observed frequency of 100 GHz and a magnetic 
field of 100 G. A critical Lorentz factor of 7crit = 10 was used to define 
the ultra-relativistic electron population. The three curves differentiate three 
assumptions for the angle between the magnetic field and the direction of 
propagation, 0. 

when — 80° and the dotted line denotes — 10°. At low 
values of the electron temperature, the magnetic field geom- 
etry is practically irrelevant with the modes remaining nearly 
purely circular 

In contrast, for a thermal plasma with a temperature kT^ > 
mc^, the modes are elliptical resulting in generalized Fara- 
day rotation. When the temperature climbs to values so that 
kT^/rac? ~ 7crit, the wave modes become linear (T = 0) and 
radiation may undergo significant conversion between linear 
and circular polarization. The magnetic field geometry has a 
large influence when the modes are elliptical. When the field 
is closely aligned to the direction of propagation, is small 
and the modes remain circular up to larger temperatures, in- 
creasing the importance of Faraday rotation. However, when 
the field is nearly perpendicular to the direction of propaga- 
tion, (f) is large, increasing the ellipticity of the wave modes 
and the importance of conversion. Both semi-analytic and 
numerical simulations of RIAFs models predict relativistic 
electron tem perature s in the inner par t of the accretion flows 
dYuanetal .1.2003: Sharmaet al.ll2()07h . Therefore, the wave 
modes will be elliptically polarized and generalized Faraday 
rotation will be important. 

3. APPLICATION TO MODELS OF Sgr A* 

Linear polarization has been observe d from Sgr A* in 
the sub-millimeter range of the spectrum (lAitken et al.ll2000l; 
iBower et all 120051; |Marrone et al ."20061). At these frequen- 
cies, the emission originates from syn chrotron emittin g ther- 
mal electrons (e.g., Ozel et al. .2000t lYuan et al.1120031) . Since 
the emission is optically thick and the synchrotron emissiv - 
ity is a very steep function of radius jLoeb & Waxmanll2007l) . 
each frequency can be mapped to a specific emission ra- 
dius in the accretion flow. Sub-millimeter VLBI observa- 
tions ha ve resolved the image of Sgr A* at a few wave- 
lengths dKrichbaum et al.li2006) . all owing a tentative calibra - 
tion of the radius-frequency relation dLoeb & Waxmanl2007h : 




101° 

Frequency (Hz) 

Fig. 3. — Axial ratio of the polarization ellipse (T; eq.[TJ at the photosphere 
for each frequency in a model of the RIAF onto Sgr A*. The three lines 
differentiate the results among three different magnetic field geometries. 

vii — 1.0rf3^, where vn = i^/(10^^ Hz) and ri3 = 
r/(10^^ cm), and a distance of 8 kpc has been assumed to 
Sgr A*. iLoeb & WaxmanI (l2007h have used this relation and 
the observed spectrum of Sgr A* to estimate the run of the 
electron temperature and magnetic field strength with radius: 



fcTe/mc^ « 14.9(2/g 



-0.4 



' 13 



and B « 21{2fgYr^^ ' G, 



where fg is a factor of order unity which we take to be 

fg = 1/ V2- These values are comparable with t hose found 
from RIAF model fits to the Sgr A* spectrum (lYuan et al.l 
l2003h . 

With these three relations, we can use the procedures de- 
scribed in § |2]to calculate the axial ratio of the normal modes 
T as a function of the radial distance from the black hole. This 
is equivalent to calculating T at the emission point for each 
frequency v. We plot in Figure [3] the results of this exercise, 
where we have calculated T from 6 to 1500 Schwarzschild 
radii, which translates to frequencies from ~ 1-100 GHz, for 
three different magnetic field geometries. We have assumed 
a black hole mass of 4 x 10^ Mg jGhez et al.ll2003l) . spheri- 
cal geometry, and have ensured that w ^ at every radius. 
The gas temperature in this model of Sgr A* is large enough 
over the entire range of radii that the plasma normal modes 
are nearly always highly elliptical and thus generalized Fara- 
day rotation will be relevent. This plot cannot give a specific 
prediction for the amount of circular and linear polarization 
expected from Sgr A* — that must wait for a full polarized 
radiative transfer calculation through an MHD accretion sim- 
ulation — but it does indicate that conversion will be impor- 
tant in the millimeter and sub-millimeter, and, thus if RIAF 
models are accurate, circular polarization is likely to be ob- 
served at these wavelengths. 

As we saw in the previous section, the magnetic field ge- 
ometry in the accretion flow will also have an effect on the 
shape of the linear modes. Figure [3] shows that a line-of-sight 
nearly aligned with the magnetic field will significantly sup- 
press conversion to circular polarization over the entire fre- 
quency range. Of course, polarization observations include 
the entire accretion flow within the beam, therefore, for a 
purely radial magnetic field, only the central region of the 
beam will observe lines-of-sight with low values of T. How- 
ever, numerical simulations of accretion flows find that the 
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velocity shear in the gas s tretches radial magnetic fields in the 
azimuthal direction (e.g- lHawley & Kro lik 2001), generating 
a toroidal geometry (i.e., a large (j>). Figure [3] shows that, in 
this case, conversion between linear and circular polarization 
will be enhanced in the beam center. It seems possible, there- 
fore, that measurements of the percentage of both linear and 
circular polarization over a range of frequency may allow a 
constraint to be placed on the overall magnetic field geometry 
in the accretion flow. 

These results clearly indicate that generalized Faraday ro- 
tation will be important in RIAFs, and thus could impact the 
interpretation of the linear polarization observations of Sgr 
A*. As discussed in §|2] if the generalized Faraday rotation is 
large enough, then the observed linear polarization is just the 
projected polarization of the plasma normal modes. Changes 
in the measured position angle with frequency are then due to 
differences in the plasma normal modes at different radii, and 
not to classical Faraday rotation. Thus, unless the generalized 
Faraday rotation is very small, the constraints on the accretion 
rate onto Sgr A* derived from the classical rotation measure 
(jMarrone et al. 2007) will not be valid. If this is the case and 
the linear polarization measurements track the changes in the 
plasma normal modes, then these observations remain power- 
ful probes of the accretion flow. 

The importance of generalized Faraday rotation in the ac- 
cretion flow onto Sgr A* can be estimated by calculating 
the angle throu gh which the polarization vector will rotate 
(lMelrosdll997h 
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where we have assumed a purely relativistic gas. Taking rep- 
resentative values and assuming ds w r we obtain 
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Therefore, it seems plausible that generalized Faraday rota- 
tion will be important in Sgr A*, and the observed linear po- 
larization measurements must be carefully interpreted. More 
precise calculations of Atp will be performed in future work. 

4. CONCLUSIONS 

The latest RIAF models of Sgr A* indicate that the thermal 
plasma is relativistic over a l arge range of radii (iYuan et al.l 
12001 iLoeb & Waxmanll2()07h . This paper has shown that, in 
general, RIAFs have highly elliptical wave modes and there- 
fore generalized Faraday rotation will dominate at millime- 
ter and sub-millimeter wavelengths. Since generalized Fara- 
day rotation in RIAFs will not depolarize the linear polarized 
component of the radiation even if the number of rotations 
is very large, linear polarization observations of Sgr A* and 
other low-luminosity accretion flows may directly probe the 
plasma normal modes, an important diagnostic of the accre- 
tion physics. In addition, the elliptical normal modes of the 
plasma will cause a cyclic conversion between linear and cir- 
cular polarization in this frequency range. As the ellipticity 
of the normal modes depend on the magnetic field geometry, 
circular and linear polarization observations of Sgr A* over a 
wide range of frequency may differentiate between radial and 
toroidal magnetic field geometry in the accretion flow. Exact 
predictions of these effects require a detailed numerical sim- 
ulation which is planned for future work. It is hoped that the 
general considerations presented here will help motivate fu- 
ture linear and circular polarization observational campaigns 
on Sgr A*. 
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